Rheology of gelling polymers in the Zimm model 
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In order to study rheological properties of gelling systems in dilute solution, we investigate the viscosity 
and the normal stresses in the Zimm model for randomly crosslinked monomers. The distribution of cluster 
topologies and sizes is assumed to be given either by Erdos-Renyi random graphs or three-dimensional bond 
percolation. Within this model the critical behaviour of the viscosity and of the first normal stress coefficient 
is determined by the power-law scaling of their averages over clusters of a given size n with n. We investigate 
these Mark-Houwink like scaling relations numerically and conclude that the scaling exponents are independent 
of the hydrodynamic interaction strength. The numerically determined exponents agree well with experimental 
data for branched polymers. However, we show that this traditional model of polymer physics is not able to 
yield a critical divergence at the gel point of the viscosity for a polydisperse dilute solution of gelation clusters. 
A generally accepted scaling relation for the Zimm exponent of the viscosity is thereby disproved. 



I. INTRODUCTION 

The influence of hydrodynamic interactions on critical 
rheological properties of gelling polymeric systems has been 
discussed controversely for many decades. In particular, the 
experimental values for the exponent fc, which governs the 
divergence of the shear viscosity ?/ cx with vanishing dis- 
tance e to the critical gel point, scatter considerably, see Ta- 
bleUfor some examples. In order to interpret the wide scatter 
of the data, they are often related to either Zimm or Rouse dy- 
namics, depending on whether hydrodynamic interactions are 
believed to be relevant or not. In this paper we intend to elab- 
orate on the validity of this interpretation, so let us be precise 
with the labels. By definition, the Rouse model' neglects both 
hydrodynamic and excluded- volume interactions. Its straight- 
forward generalization from Unear polymers to a gelling melt 
of randomly crosslinked monomers provides a microscopic 
framework, within which one can derive an exact scaling re- 
lation for the viscosity exponenti2iii^ The Zimm model,^ by 
definition, takes hydrodynamic interactions into account on a 
preaveraged level, but still neglects excluded-volume interac- 
tions. We are not aware of a microscopic approach based on 
the Zimm model which allows for an exact analytic compu- 
tation of the viscosity exponent for a gelling polymeric solu- 
tion. Other models for gelling polymers, which go beyond 
the Rouse or Zimm model by incorporating excluded-volume 
effects or fluctuating hydrodynamic interactions, are conjec- 
tured to belong to different universality classes and will not be 
considered here. 

Scaling theorjii^iiii^ has proven to be a powerful tool to 
describe the properties of polymeric systems. The relaxation 
time tn of a typical cluster of n monomers (henceforth n will 
be referred to as the size of the cluster) is estimated to be 
tn ^ R^/Dn, where i?„ is the radius of gyration and £)„ 
the diffusion constant of the cluster The scaling i?„ ~ n^^'^f 
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of the radius of gyration of a cluster with size n is determined 
by the Hausdorff fractal dimension d / . The diffusion constant 
is assumed to scale like Z?„ ^1/ i?„ in the Zimm model. This 
assumption is based on the Stokes Einstein relation which is 
valid for linear polymers and still holds for the diffusion of 
fractal polymer clusters in the Zimm model.''' The average 
contribution of clusters of size n to the viscosity is then given 
by ?7„ ^ tn/n. This implies the scaling 
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d/df - 1 



(1) 



TABLE I: Experimental data for the critical exponent k of the vis- 
cosity at the gelation transition. 



where d is the spatial dimension. With an underlying distribu- 
tion of cluster sizes, which is widely believed to follow the 
scaling laws of percolation, this gives rise to the exponent 
k = {1 — T + d/df)/a for the averaged viscosity in terms 
of the static percolation exponents. 

Beside critical properties, recent publications aim at the 
dynamics of single clusters with particular topologies within 
the Zimm model. Refs. and |3 consider the Zimm d^ 
namics of star-shaped clusters and dendrimers, and Ref. (fj 
analyses the relaxation behaviour of fractal (Sierpinski-type) 
clusters in the Zimm model. The latter authors mention the 
possibility of non-universal behaviour The question of non- 
universality is also raised from computer simulations^-' of 
gelling liquids under the influence of solvent particles. 

In this paper we investigate the viscosity and the nor- 
mal stresses in the Zimm model for randomly crosslinked 
monomers. The distribution of cluster topologies and sizes is 
assumed to be given either by Erdos-Renyi random graphs or 
three-dimensional bond percolation. The details of the model 
are described in SectionUl] Within this model, the critical be- 
haviour of the viscosity 77 and of the first normal stress coef- 
ficient \I'('^^ is determined by the scaling with n of the partial 
averages ?]„ ^ n'''', respectively "ifiP ^ n''*, over clusters 
of size n (Section HIH. In Section [TV] we investigate these 
Mark-Houwink like scaling relations numerically for differ- 
ent strengths of the hydrodynamic coupling constant. We con- 
clude in Section|V]that (i) these scaling relations are governed 
by universal exponents 6,, and 65,. This conclusion is substan- 
tiated by comparing our results to those for ring polymers in 
the Zimm model, which are known to exhibit universal be- 
haviour, (ii) We find that the scaling relation Q does not 
agree with our numerical data and, hence, does not describe 
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the viscosity in the Zimm model for randomly crosslinked 
monomers. 



dynamics is recovered, if the non-diagonal terms i ^ j of the 
mobility matrix are neglected. The Gaussian thermal-noise 
force fields f, (r, t) in (|3 have zero mean and covariance 



II. DYNAMIC MODEL AND ITS SOLUTION 



f,(r, t) f J (r', t') = 2 H,,, (r - r') Sit - t') 



(6) 



A. Hydrodynamic Interactions 

We consider N point-like monomers, which are charac- 
terized by their time-dependent position vectors Ri(t), i = 
1, . . . , A^, in three-dimensional Euclidean space. Permanently 
formed crosslinks constrain M randomly chosen pairs of par- 
ticles {ie,je), e = 1,...,M. We study the dynamics of 
crosslinked monomers in the presence of a solvent fluid, giv- 
ing rise to hydrodynamic interactions between the monomers. 
Purely relaxational dynamics in an incompressible fluid sub- 
jected to an external space- and time-dependent flow v(r, t) is 
described by the equation of motionSiS^ 



Here S stands for the Dirac-delta function and the overbar in- 
dicates the Gaussian average over all realizations of f . 

In order to determine the model completely, it only re- 
mains to specify the probability distribution of the crosslink 
configurations. We shall discuss two different types of proba- 
bility distributions: (i) crosslinks are chosen independently 
with equal probability for every pair of monomers, corre- 
sponding to Erdos-Renyi random graphs,^** and (ii) a distri- 
bution of crosslinks, which generates clusters amenable to the 
scaling description of finite-dimensional percolationJS. The 
precise characterization of these distributions is given below. 
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(2) 



Here, crosslinks are modelled by Hookean springs in the po- 
tential energy 

g Af 2 3^ 



(3) 



where the length a > plays the role of an inverse crosslink 
strength and physical units have been chosen such that U-qT — 
1. A given crosslink configuration Q = {ie,je}eLi is speci- 
fied by its iV X A^-connectivity matrix T. Moreover, we impose 
a simple shear flow 



v(r,i) 




(4) 



which is characterized by its time-dependent shear rate j{t). 
The mobility matrix is given by 



H.,,(r) 
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(5) 



The diagonal term in (jSj accounts for a frictional force with 
friction constant ( that acts when a monomer moves relative 
to the externally imposed flow field 0. The non-diagonal 
term reflects the influence of the motion of monomer j on 
the solvent at the position of monomer i and is given by the 
Oseen tensor Here rjs denotes the solvent viscosity, Sij 
the Kjonecker symbol, 1 the three-dimensional unit matrix 
and the dagger indicates the transposition of a vector Rouse 



B. Preaveraging Approximation 

The equation of motion (|2ji is nonlinear due to the nonlin- 
ear dependence of the mobility on the particles' positions. A 
simple but uncontrolled approximation is the so-called preav- 
eraging approximation that was first introduced by Kirkwood 
and Riseman^^ and Zimm.^ In this approximation the mobility 
matrix (|5} is replaced by its expectation value (Hi j)cq, which 
is computed with respect to the equilibrium distribution, i.e. 
the Boltzmann weight ~ e^^. Due to rotational invariance of 
the potential (|3}, the averaged mobility matrix is a multiple of 
the identity matrix (Hij (Ri — Rj))cq = H""! 1, where 



R,; — R,- 



(7) 



cq 



In the computation of 0, care has to be taken of the zero 
eigenvalues of the connectivity matrix, corresponding to the 
translation of whole clusters. To this end we regularize the 
potential Q by adding a confining term 3uj/ (2a^) X^iLi ' 
Ri and letting ui > tend to zero subsequently. The average 
in (0 is conveniently performed via the Fourier representation 
of 1/ 1 r I , and the result 



lim ( [GHIm + [GH],,, 



2[G(^)] 



-1/2 



(8) 



involves the resolvent G{u!) :— (F + ujI)^^ of F. The limit 
i is taken by expanding the resolvent G(aj) = Eq/u + 
Z + 0{uj) in terms of uj. Here Z := (1 - Eq) /F is the Moore- 
Penrose inverse of the connectivity matrix, i.e. the inverse of 
F restricted to the subspace of non-zero eigenvalues. More- 
over, 1 denotes the N x A^-unit matrix and Eq the projector 
on the nullspace of F, which is spanned by the vectors that 
are constant when restricted to any one cluster of crosslinked 
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monomers. More precisely, the matrix element [Eo]i.j is given 
by the inverse number of monomers of the cluster if i and j 
are in the same cluster and zero otherwise (cf. Sec. II. D in 
Ref.|3|for details). Hence, the right-hand side of (|8j vanishes 
for w I whenever i and j belong to different clusters. Con- 
sequently, the preaveraged mobility matrix shows corre- 
lations of different particles only if these particles are in the 
same cluster, in other words it is block-diagonal and within 
one block given by 



^,3 



5ij + (1 - 5i,j)h [k^ T^lT^i,j) 



(9) 



For convenience we introduced the function h{x) = -^/x/tt 
and the quantity TZi^j := Z^.i + j — 2Zi j, which can 
be interpreted as the resistance between nodes i and j in 
a corresponding electrical resistor network.^' The parameter 
K :— y/WpK C,/ {Q-nrisO) plays the role of the coupling con- 
stant of the hydrodynamic interaction. Note that this defi- 
nition of K differs from that of other authors by a factor of 
\/2,^^ respectively VG/tt.^^ Formally setting k = in (|9j 
yields V\f^- = (^^Sij, and the Zimm model for gelation 
reduces to the Rouse model for gelationi ^'-^'^^i^i^i^-^ It is well 
known that the Oseen tensor does not give rise to a positive- 
definite mobility matrix for all possible spatial configurations 
of monomers. This defect is cured if the Rotne-Prager- 
Yamakawa tensos^i^is used instead. Again, the preaveraging 
procedure is done with a confining potential which is switched 
off afterwards. The function h is then given by"*^ 



h{x) = erf{y/x) 



1 1 — exp(— x) 



(10) 



It involves the error function erf (x) and recovers the form of 
the preaveraged Oseen-Tensor asymptotically as a; | 0. As 
a result of preaveraging we obtain the Zimm model for cross- 
linked monomers in solution 



The resulting equation of motion for Rj {t) coincides with that 
of the Rouse model for crosslinked monomers, ^i-^i^'^'^'-^-^ if one 
replaces the connectivity matrix T by 



/2 



(14) 



in the latter Different coordinate transformations are com- 
monly used to establish this formal relation between the two 
models. We prefer (I13> . because then the transformed equa- 
tion of motion involves the symmetric matrix (I14> . The result- 
ing monomer trajectories for (transformed) initial data [to) 
are therefore given by 



to 



(15) 



as follows e.g. from Sec. II. C in Ref. 3. The solution il5\ is 
expressed in terms of the transformed thermal noise with zero 
mean and covariance 



(16) 



and the time evolution in the simple shear flow Q is charac- 
terized by the N x A^-matrix 



and the 3 x 3-matrix 



1 f^,dsjis) 







(17) 



(18) 
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-R,(t)-v(R,(t),t) = -^H 
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Here, the covariance of the thermal noise is given by 



(11) 



(12) 



Since both the connectivity matrix F and the preaveraged mo- 
bility matrix H'^'i are block-diagonal, it follows that clusters 
move independently of each other in this model. 



C. Formal Solution 

The Zimm equation (II 1> is linear, hence it can be solved 
exactly. This is most conveniently done by introducing new 
coordinates Rj (t) through the coordinate transformation 



N 

E 



[(H°^) 



1/2 



(13) 



Finally, the solution of the Zimm equation ( II W is obtained by 
inserting dlSt in ( I13t . 



III. OBSERVABLES 



A. Shear Stress 



We shall focus on the viscosity r] and the first and second 
normal stress coefficients ^'(^^ and ^f^^^, respectively. There- 
fore we need to compute the intrinsic shear stress cr(i) as a 
function of the shear rate 7(t). Following Chap. 3 in Ref.l24l 
or Chap. 16.3 in Ref.|3 we express the shear stress in terms 
of the force per unit area exerted by the polymers 



N 

lim -^^F,(<)RlW. 



(19) 



Here, Ri(t) is the solution of the equation of motion dl 1> with 
some initial condition Ri(to) at time in the distant past (so 
that all transient effects stemming from the initial condition 
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have died out). Moreover, po stands for the monomer con- 
centration and Fi{t) := ~dV/dIli{t) is the net spring force 
acting on monomer i at time t. Using the transformation il3\ 
and the solution (I15> . it is readily shown^ that the stress ten- 
sor ( I19> is given by 

<T(t) =x(o)i+ r dt' x{t~t')j{t') 

-/ — oo 

/ 2j^Usi{s) 1 



V 




(20) 



for arbitrary strengths of the shear rate 7(t). Here, we have 
defined the stress-relaxation function 



(1 - Eo) exp 



6t 



(21) 



as a trace over the subspace of non-zero eigenvalues of F. 

For a time-independent shear rate 7, the shear stress (l20l 
is also independent of time. The (intrinsic zero-shear) viscos 
ity T] is then related to shear stress via 



ipo 



and the normal stress coefficients are defined by 



7V0 



15(2) 



respectively. Hence, the viscosity \22\ is given by 



viG) = - dixit) = - — Tt 
Po Jo 3 2N 



(22) 



(23) 



(24) 



for a fixed realization Q of crosslinks. It is determined by the 
trace of the Moore-Penrose inverse of T{Q). According to 
(|20j and ( I23> . the second normal stress coefficient ^f'^^ van- 
ishes always, whereas 



vl/(i)(g) = A 

Pq Jo 



1-Eo 



(2^ \ 2 -j^ 

dttxit) = ( — I Tr _ 



(25) 



Again, we have made explicit the dependence on Q in J25t . 
This will be convenient for computing the average over all 
crosslink reaUzations in the next subsection. 



B. Disorder Average and Critical Beliaviour 

Each crosslink reaUzation Q defines a random labelled 
graph on the set of monomers, which can be decomposed into 
maximal path-wise connected components or clusters 



K 



(26) 



Here, Afk denotes the A:-th cluster with Nk monomers out of a 
total of K clusters (all depending on Q). We also refer to Nk 
as the size of the^cluster JVk- The associated modified con- 
nectivity matrix F from ( I14t is of block-diagonal form with 
respect to the clusters. Therefore one can decompose any 
observable of the type A{g) = N-^TTf{f{g)), where / is 
some function on the reals, into contributions from different 
clusters according to 



K 



fe=l 



iV 



(27) 



Here, we have defined A{J\fk) := jVr^Tr/(r(7Vfc)). In par- 
ticular, J27> holds for the viscosity J24li and for the first normal 
stress coefficient (l25t . 

In order to compute the average (A) of the observable 
A over all crosslink realizations in the macroscopic limit 
M 00, ^ 00 with fixed crosslink concentration 
c := M/N, we have to specify the statistical ensemble that 
determines the realizations of crosslinks. Two distributions 
of crosslinks will be considered, (i) Each pair of monomers 
is chosen independently with equal probability c/N, corre- 
sponding to Erdos-Renyi random graphs, which are known 
to resemble the critical properties of mean-field percolation. 
After performing the macroscopic limit, there is no macro- 
scopic cluster for c < Ccrit = 1/2 and almost all clusters 
are trees. Furthermore, all n"^^ trees of a given size n are 
equally likely, (ii) Clusters are generated according to three- 
dimensional continuum percolation, which is closely related 
to the intuitive picture of gelation, where monomers are more 
likely to be crosslinked when they are close to each other. 
Since continuum percolation and lattice percolation are be- 
lieved to be in the same universality clasSjS we employ the 
scaling description of the latter. It predicts2S a cluster-size 
distribution of the form 



K 



k=l 



^ exp{— 7i/n*} 



(28) 



for e := (ccrit — c) ^ 1 and n 00 with a typical cluster 
size n*{e) ^ e^^^°' that diverges as e ^ 0. Here, a and r are 
(static) critical exponents. 

For the computation of the average (A) over all crosslink 
realizations Q it is convenient to introduce partial averages 



fc=i 



(29) 



of A over all clusters of a given size n. Using (|^} and re- 
ordering the clusters, one gets the identity 



K 



J2^AiMk))=J2nr„{A)„, 



k=l 



N 



(30) 



which is valid in the absence of an infinite cluster Now sup- 
pose one has the scaling 



fc=i 



Ar. 



e=0 



(31) 
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of the partial average at the critical point as n ^ cjo. Due to 
the absence of relevant scales at the critical point, this is quite 
a natural behaviour. Then, (I30> and ( I28> imply the critical 
divergence 

(.4)-e-" as eiO, with u^{2-T + h)/(j (32) 

for the crosslink-averaged observable A, provided that m > 0. 
We will therefore study the scaling 

r]n ^ n''" and - n''* (33) 

as 71 ^ oo to explore critical rheological behaviour at the 
gelation transition within the Zimm model. 

Formulas like ( I30> - ( I33> may also be familiar from scal- 
ing theories for gelation. We go beyond such approaches in 
that we have mapped the dynamical properties of a gelling 
molecular system to a percolation problem, see e.g. MA\ and 
M5\ . This mapping has been fully derived within a (semi-) 
microscopic dynamical model, the Zimm model for randomly 
crosslinked monomers, and not merely postulated from ad- 
hoc assumptions, as is usually done in scaling theories. In the 
following section we describe the numerical solution of the 
percolation problem. 



IV. NUMERICAL RESULTS 
A. Erdos-Renyi Random Graphs 

For numerical purposes it is convenient to compute the 
eigenvalues of the non-symmetric matrix F := V\'^'^V rather 
than those of f (H<=q)i/2r(H<=q)i/2 because this prevents 
us from computing the expensive square root of H'^'^. The fact 
that F and F have the same eigenvalues can easily be proven 
by observing that if is an eigenvector of F with correspond- 
mg eig envalue A then (H'=q)±i/2 i/' is a right/left eigenvector 
of F with the same eigenvalue A. 

As already mentioned in Sec. 1111 Bl the average (•)„ ex- 
tends over all n"^'^ equally weighted labelled trees of size n 
in the case of Erdos-Renyi random graphs and, hence, is in- 
dependent of the crosslink concentration c. Random labelled 
trees of a given size have been generated via the Priifer algo- 
rithm and handled with the LED A library^ The preaveraged 
mobility matrix (|9} is computed with the function h from ( II OK 
corresponding to the Rotne-Prager-Yamakawa tensor The 
resistances TZij in trees reduce to shortest paths, that is graph 
distances, which are calculated with the Dijkstra algorithm.^** 
The eigenvalues of F are then computed with the LAPACK 
library. For suitable, logarithmically equidistant cluster sizes 
71 G [2, 4000] we average the viscosity and the normal stress 
coefficient over 50 trees, which turned out to yield an accept- 
able computer-time/accuracy trade-off. In Figs. [Qa) and (b) 
we plot rjn and "in^ as a function of ti on a double-logarithmic 
scale for different values of the hydrodynamic interaction pa- 
rameter K,. According to ( I33l l the exponents 6,, and are 
obtained by power law fits in the large 77-range, for which we 



choose the interval n S [700, 4000], see Fig.^c). For the vis- 
cosity the exponent decreases from 6,, = 0.28 for k = 0.05 
to 6,, = 0.11 for K ~ 0.3. The Rouse exponent for k ~ 
is exactly given by"* 6,, = 1/2. The exponent 6^ of the nor- 
mal stress coefficient ranges from 6* = 1.2 for k = 0.05 to 
~ 0.73 for K — 0.25. The Rouse value for k = is exactly 
given"* hy b^^ ~ 2. 



B. Three-Dimensional Percolation 

For the generation of clusters according to three- 
dimensional bond percolation we apply the Leath 
Algorithm.-''' It generates a sequence {JVi}f"^i of clus- 
ters, in terms of which the disorder average is readily 
computed via (A) = limL^oo ^^^^ A{J\fi). This 

implies {A),, = limi^oo E/=i SN,,nA{M)/ ^N,,n for 
the average over clusters of size n. The algorithm has been 
tested by verifying the scaling of the cluster-size distribution 
T„. Second, for small values of n, we compared the number 
of clusters with known exact values4S Third, we verified that 
the exponent 2/df, which governs the scaling of the squared 
radius of gyration as a function of cluster size n at the critical 
pointp^ comes out as 2/2.53 from the simulation. For each 
generated cluster the resistances TZij are computed from the 
Moore-Penrose inverse Z of the connectivity matrix F - see 
below Eq. (|9} - and inserted into (|9} with h corresponding 
to the Rotne-Prager-Yamakawa tensor. The eigenvalues of 
r are then computed with the LAPACK library. We were 
forced to restrict cluster sizes to values n < 4000 due to the 
limited amount of memory, which is required for the gener- 
ation and diagonalization of the matrix product F = H'^'^F. 
Moreover, for calculating disorder averages we restrict the 
number of realizations pertaining to a given cluster size to 
a maximum of 50. However, within the present numerical 
effort this maximum number is not even attained for larger 
cluster sizes. Therefore the disorder averaged quantities are 
still subject to fluctuations. In order to obtain smooth curves 
for rjn and we have also smoothed out the raw data by 
performing a running average over cluster sizes in the window 
[n — 5,n + 5]. The thus obtained values for rjn and 'i'n are 
plotted in Figs.QJd) and (e), respectively, as a function of n on 
a double-logarithmic scale for different values of k. The ex- 
ponents b,j and 6^, extracted by fitting the curves in Figs.[3d) 
and (e) to a power law in the interval n e [800,4000], are 
shown in Fig. [flf)- The numerical values for 6^ are nearly 
identical to those obtained for Erdos-Renyi random graphs. 
Again, one observes a decrease from 6^ — 0.21 for k = 0.05 
to 5^ = 0.11 for K — 0.3. The exponent 6* of the normal 
stress coefficient ranges from b^ = 1.1 for k — 0.05 to 
6* = 0.78 for K — 0.25. The corresponding Rouse values 
for K = follow from exact analytical arguments^ii^ and are 
given by 6,, = (2/4) - 1 « 1/2 and 6v[, = (4/ds) - 1 « 2, 
respectively. Here, ds « 4/3 is the spectral dimension of the 
incipient percolating cluster, whose numerical value is very 
well approximated by the Alexander-Orbach conjecture ^22. 
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FIG. 1 : Numerical data to determine the scaling <33t for random clusters in the case of Erdos-Renyi random graphs (left column) and three- 
dimensional bond percolation (right column). In each case the averaged viscosity r]„ (top) and normal stress coefficient ^'Ji^' (middle) are 
plotted for different strengths of the hydrodynamic interaction parameter «; as a function of the cluster size n on a double logarithmic scale. 
Power-law fits to the data yield the exponents 6,, and 6* as a function of n (bottom). 



C. Ring Polymers 

We suspect that the observed variation of the exponent 
values with k may be due to crossover and finite-size ef- 
fects. To clarify this question it is useful to study a system 



where the exponents are known analytically. Therefore we 
(re-)investigate the viscosity Turing and the first normal stress 

coefficient ^'^iilg of ring polymers in the Zimm model with the 
Rotne-Prager-Yamakawa tensor The scaling of both quanti- 
ties with ring size n as n ^ oo can be deduced from long- 
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standing analytical results, - which lead to b 



1/2 and 



6* : 



2. We focus here on the onset of this asymptotic 



behaviour and how it is affected by crossovers for different 
K. This provides us with a reference system when discussing 



the scaling of ?7„ and "i! 
Section|Vl 



(1) 



in the case of random clusters in 



Due to the cyclic structure of a ring polymer the associ- 
ated matrices H'^'^ and F are circulant matrices. Hence, they 
are simultaneously diagonalizable. In fact, the j-th compo- 
nent of the l-th eigenvector of F for a ring of size n is ex- 
plicitely given by = exp(i27rj7/ri), and as a result the 
eigenvalues can be written in terms of Fourier transformsi^ 
Therefore, r/^ing and '^l^lg are efficiently computed by Fast 
Fourier Transformation up to ring sizes n — 10^. The result- 
ing viscosity and the first normal stress coefficient are shown 
m Figs.lSa) and (b) on a double logarithmic scale. The data 
is then fitted to a power law in two different fit ranges. In ad- 
dition to a fit in the terminal large-n range, n G [10'', 10^], we 
performed a second fit in the range n S [500, 5000], which is 
roughly where we had to do the fits in the random-cluster case. 
The fit exponents are shown in Fig.|3c). Apparently, they de- 
pend on the fit range. For k — 0.05 we find fe,,,ring = 0.69 
from the small-n fit. This value clearly exceeds the theo- 
retical one 6,,, ring — 1/2. Even the corresponding value 
^.ring = 0.58 from the large-n fit still has an error of 36%. In 
contrast, for k — 0.3 both values, fe,),ring — 0.51 and 0.50, are 
quite close to the exact one. 

In fact, given the Fourier representation of the eigenvalues 
of F, it is straightforward to demonstrate the occurrence of a 
crossover at n « tt/k^ from Rouse behaviour. Typing ^ n, 
to the asymptotic Zimm behaviour 77ring ~ v}!'^ j k for all 
n ^ kt"^ . Hence, the larger k, the less important is resid- 
ual Rouse behaviour in the numerical data for the scaling of 
77iing. The same holds true for 4'|j'j|g. Unfortunately, choosing 
larger values for n is not a practicable way out for obtaining 
good-quahty data. This is because for large k the asymptotics 
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of ( I10> as a; — > oo becomes noticible and leads to the transient 
behaviour raring ~ kvl-' for intermediate n. We have observed 
such a behaviour for (unphysically large) k > 10 (not shown). 
But even the data for k = 0.5 and k = 1.0 in Fig. |2] are still 
slightly influenced by ( I34> . 

In summary, whereas there is a generic overestimate of 
the scaling exponents for small k due to residual Rouse be- 
haviour, the exponents are underestimated for higher k due 
to the asymptotics ( I34l i. The optimal value for minimal finite- 
size effects in ring polymers appears to be k w 0.3 in Fig.|2jc). 



FIG. 2: Numerical data to determine the scaling <33t for ring poly- 
mers. The viscosity raring (a) and the normal stress coefficient 'I'Jif^g 
(b) are plotted for different strengths of the hydrodynamic interac- 
tion parameter k as a function of the cluster size n on a double 
logarithmic scale, (c) shows the exponents 6,7, ring and 6*, ring from 
power-law fits to the data of (a) and (b). The fits were performed 
for two different ranges of cluster sizes n. Additional data points for 
K = 0.05, 0.15 and 0.25 in (c) stem from curves which have been 
omitted in (a) and (b) for reasons of clarity. The two horizontal lines 
indicate the exact values br),ring = 1/2 and 6*, ring = 2. 
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V. DISCUSSION 

Using the Zimm model for randomly crosslinked 
monomers, we have determined the scaling ( I33> of the av- 
eraged viscosity and of the averaged first normal stress coeffi- 
cient over clusters of a given size n. Figs.fHc) and (f) display a 
crossover from the Rouse values at k = to the Zimm values 
at non-zero k. We estimate the latter as 

fe,, wO.ll and b^fiiO.77, (35) 

from our data for k = 0.3. A detailed discussion of this choice 
of K and of possible origins of the dependence of bjj and 6* 
on K will be given below. Within the accuracy of our data, the 
exponents are the same for both Erdos-Renyi random graphs 
and three-dimensional percolation. 

The critical behaviour of the averaged viscosity (77) ^ 
and of the averaged first normal stress coefficient 
^ for a polydisperse gelling solution of cross- 
linked monomers then follows from (I35> and (I32> . For the 
viscosity this implies a finite value at the gel point for both, 
Erdos-Renyi random graphs and three-dimensional bond per- 
colation. In contrast, the first normal stress coefficient is 
found to diverge with an exponent that depends on the cluster 
statistics. Choosing the cluster statistics according to Erdos- 
Renyi random graphs, we find £ w 0.54. The case of three- 
dimensional bond percolation leads to the higher value I w 
1.3. These exponent values are less than a third in magnitude 
than the corresponding exact analytical predictions £ = 3, re- 
spectively £ K, 4.1 of the Rouse model for randomly cross- 
linked monomers"^-^ '' with the corresponding cluster statistics. 

The dependence of the critical exponents b-q and fe* on 
the hydrodynamic interaction strength k in Figs.^^c) and (f) 
may be due to finite-size effects. In particular the onset of 
the true asymptotic regime of these quantities may depend on 
K. In order to better understand finite-size effects, we have 
examined the Zimm dynamics of polymer rings in Sec. lIV Cl 
and determined the scaling of the viscosity ?7ring ~ n^i-^'^s 
and of the first normal stress coefficient ^^lUg ^ ^^''f ring ^yjtjj 
the ring size n. For rings one can access much higher values 
of n as for random clusters, see Figs.|2ja) and (b). In partic- 
ular, the exactly known scaling exponents fo,,,ring = 1/2 and 
&*,ring — 2, which are universal in k > 0, can be extracted 
from our data in Fig.lJJc). However, if we did not exploit the 
full range of available ring sizes and restricted the fit to those 
lower values of n which could also be accessed for random 
clusters, then universality would be veiled by finite-size ef- 
fects. Finite-size effects are more pronounced for k < 0.15 
and K > 0.5. Thus, we conclude (i) that the random-cluster 
data have not reached either the asymptotic large-ri regime yet 
for K < 0.15 in Fig.^ (ii) that the asymptotic regime is uni- 
versal and (iii) that the data for k — 0.3 should be the most 
reliable ones. 

The exponent 6,, has also been investigated experimen- 
tally. In Ref. ^ measurements on randomly branched 
polystyrenes have been performed, resulting in 6^ e 



[0.2, 0.25]. Measurements on branched polyethyleneimine — 
yield the slightly higher value w 0.31. Brownian-dynamics 
simulations of hyperbranched polymers were performed in 
Ref. 44. They also account for fluctuating hydrodynamic in- 
teractions corresponding to k = 0.35, as well as for excluded- 
volume interactions and lead to 5^ = 0.13. This result is re- 
markably close to our finding 6^ w 0.11 for the highest cou- 
pling strength k — 0.3 that we have considered, whereas the 
experimental findings are consistently above our value (see 
the discussion below). 

Next we compare our findings with the scaling argu- 
ment which is summarized in Eq. Q. For phantom clus- 
ters, i.e. in the absence of excluded-volume interactions, the 
Hausdorff fractal dimension is equal to the Gaussian fractal 
dimensioni^il^ d^^'' :— 2ds/{2 — ds), where dg is the spec- 
tral dimension. Here we estimate ds ~ 4/3 according to 
the Alexander-Orbach conjecture, which is known to be an 
excellent approximation albeit not being exact. For 0? = 3 
Eq. Q then implies 6^ « —1/4, and for d — 6 one has 
b,f « 1/2. The latter value corresponds to Erdos-Renyi ran- 
dom graphs, whose critical properties are identical to those of 
mean-field percolation.^^ Both values can be definitely ruled 
out by our data. Thus we conclude that the scaling relation 
Q does not apply to the Zimm model for randomly cross- 
linked monomers. This failure comes as a surprise because 
it is known from a recent investigation of diffusion constants 
within this model that the exact results are in accordance 

(G) 

with long standing scaling relations when inserting dK ' for 
df. 

Coming back to the experimental fc-values listed in Ta- 
ble m and considering also the exact prediction k = {1 — 
T + 2/ds)/a w 0.71 of the Rouse model for gelling 
monomers, ^'^'^ we conclude that an explanation for the broad 
scatter of the data in the literature calls for additional relevant 
interactions than those accounted for in the Zimm or Rouse 
model. This may be due to the preaveraging approxima- 
tion. In particular, it throws away hydrodynamic interactions 
among different clusters. But we do not expect this to be the 
sole relevant simplification of the Zimm model, because linear 
polymers show a decrease in the viscosity when abandoning 
the preaveraging approximation,^^ and effects of preaveraging 
for branched molecules are even more pronounced than those 
for linear ones4i Rather it seems that there are no satisfac- 
tory explanations without considering excluded-volume inter- 
actions. Indeed, simulations'*** of the bond-fluctuation model 
deliver higher values few 1.3 in accordance with the scaling 
relation k = 2v ~ which arises from heuristically merging 
Rouse-type and excluded-volume properties. 
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